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A class of exact solutions of the gravitational field equations in the vacuum on the brane are 
obtained by assuming the existence of a conformal Killing vector field, with non-static and non¬ 
central symmetry. In this case the general solution of the field equations can be obtained in a 
parametric form in terms of the Bessel functions. The behavior of the basic physical parameters 
describing the non-local effects generated by the gravitational field of the bulk (dark radiation and 
dark pressure) is also considered in detail, and the equation of state satisfied at infinity by these 
quantities is derived. As a physical application of the obtained solutions we consider the behavior 
of the angular velocity of a test particle moving in a stable circular orbit. The tangential velocity 
of the particle is a monotonically increasing function of the radial distance and, in the limit of large 
values of the radial coordinate, tends to a constant value, which is independent on the parameters 
describing the model. Therefore a brane geometry admitting a one-parameter group of conformal 
motions may provide an explanation for the dynamics of the neutral hydrogen clouds at large 
distances from the galactic center, which is usually explained by postulating the existence of the 
dark matter. 
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I. INTRODUCTION 

In two recent papers Q. (I several classes of conformally symmetric solutions of the static gravitational field 
equations in the brane world scenario (S'], in which our Universe is identified to a domain wall in a 5-dimensional anti- 
de Sitter space-time, have been obtained. The solutions have been derived by using some methods from Lie group 
theory. As a group of admissible transformations the one-parameter group of conformal motions has been considered. 
More exactly, as a starting point it has been assumed that the metric tensor g^ on the brane has the property 
Lcg^ = ip (r) g^v , where the left-hand side is the Lie derivative of the metric tensor, describing the gravitational field 
in vacuum on the brane, with respect to the vector field and %p, the conformal factor, is an arbitrary function of 
the radial coordinate r. 

The assumption that in the brane world model the vacuum outside a matter distribution has a self-similar structure 
is not at all arbitrary; it has a strong theoretical and even observational basis. The assumption that the space-time 
admits, beside the spherical symmetry, a one-parameter group of conformal motions, is a particular case of the 
geometric self-similarity, a property of the metric which has extensively been investigated in the literature (for a 
review of this important subject see [|j). Geometric similarity should be distinguished from physical similarity, which 
is a property of the matter fields. These two properties are not necessarily equivalent. A ’’similarity hypothesis”, 
which assumes that in a variety of physical situations, solutions of strongly non-linear differential equations may 
evolve naturally to a self-similar form, even if they start out in a more complicated form, has also been proposed |5| . 
The expansion of the Universe from the big bang and the collapse of a star to a singularity tend to self-similarity 
in some circumstances. Spherically symmetric cosmological fluctuations might naturally evolve from complex initial 
conditions via the Einstein equations to a self-similar form. Self-similar asymptotics can be obtained in a wide range 
of contexts in fluid dynamics and for orthogonally transitive Gi cosmological models. The presence of the voids in 
universe and the large scale structure can be described by self-similar asymptotic solutions of the Friedmann equations 
pj. All hypersurface homogeneous locally rotationally symmetric spacetimes, which admit conformal symmetries 
have been determined in [7|, and the symmetry vectors were explicitly obtained. Hence conformal symmetry is 
generically possible, it is a mathematical property of the solutions of many non-linear differential equations and in 
some pliysical/astrophysical phenomena this property has been effectively observed. 
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Due to the correction terms coming from the extra dimensions, at very high energies significant deviations from the 
Einstein theory occur in the brane world models [g] (for a recent review of the brane world theories see 0)- Gravity 
is largely modified at the electro-weak scale 1 TeV. The cosmological implications of the brane world theories have 
been extensively investigated in the physical literature. 

For standard general relativistic spherical compact objects the exterior space-time is described by the Schwarzschild 
metric. In the five dimensional brane world models, the high energy corrections to the energy density, together with 
the Weyl stresses from bulk gravitons, imply that on the brane the exterior metric of a static star is no longer the 
Schwarzschild metric M, lll| . The presence of the Weyl stresses also means that the matching conditions do not have 
a unique solution on the brane; the knowledge of the five-dimensional Weyl tensor is needed as a minimum condition 
for uniqueness. The static vacuum gravitational field equations on the brane depend on the generally unknown Weyl 
stresses, which can be expressed in terms of two functions, called the dark radiation U and the dark pressure P terms 
(the projections of the Weyl curvature of the bulk, generating non-local brane stresses) 0, fg]. 

Generally, the vacuum field equations on the brane can be reduced to a system of two ordinary differential equations, 
which describe all the geometric properties of the vacuum as functions of the dark pressure and dark radiation terms 
01. In order to close the system of vacuum field equations on the brane a functional relation between these two 
quantities is necessary. Hence a first possible approach to the study of the vacuum brane consists in adopting an 
explicit equation of state for the dark pressure as a function of the dark radiation. This method has been systematically 
considered in 0, where some classes of exact solutions of the vacuum gravitational field equations on the brane have 
been derived and a homology theorem for spherically symmetric vacuum on the brane has also been proven. 

By using various assumptions and ansatze several classes of exact vacuum or interior solutions of the gravitational 
field equations on the brane have been obtained in 0] (the solution has the mathematical form of the Reissner- 
Nordstrom solution, in which a tidal Weyl parameter plays the role of the electric charge of the general relativistic 
solution), in (the exterior solution also matches a constant density interior) and in Il2j . 

A different approach, which avoids the consideration of some ad hoc equations of state for the dark pressure consists 
in assuming that the vacuum brane has some particular symmetries, and to investigate these symmetries by using 
Lie group techniques 0, 0. As for the vector field generating the symmetry one can assume that it is static and 
spherically symmetric 0 or non-static and spherically symmetric [2|. In both cases the gravitational field equations, 
describing the static vacuum brane with geometric self-similarity can be integrated in Schwarzschild coordinates, 
and several classes of exact solutions, corresponding to a brane admitting a one-parameter group of motions, can be 
obtained. The general solution of the field equations depends on some arbitrary integration constants. The main 
advantage of imposing geometric self-similarity via a group of conformal motions is that this condition also uniquely 
fixes the mathematical form of the dark radiation and dark pressure terms, respectively, which describe the non-local 
effects induced by the gravitational field of the bulk. Thus there is no need to impose an arbitrary relation between 
the dark radiation and the dark pressure. 

As a possible physical application of the conformally symmetric solutions with non-static symmetry of the spherically 
symmetric gravitational field equations in the vacuum on the brane the behavior of the angular velocity Vt g of the 
test particles in stable circular orbits has been considered in 0. In this case the conformal factor ij), together with 
two constants of integration, uniquely determines the rotational velocity of the particle. In the limit of large radial 
distances and for a particular set of values of the integration constants the angular velocity tends to a constant 
value. This behavior is typical for massive particles (hydrogen clouds) outside galaxies 0, and is usually explained 
by postulating the existence of the dark matter. Thus, the rotational galactic curves can be naturally explained in 
brane world models without introducing any additional hypothesis. The galaxy is embedded in a modified, spherically 
symmetric geometry, generated by the non-zero contribution of the Weyl tensor from the bulk. The extra-terms, which 
can be described in terms of the dark radiation term U and the dark pressure term P, act as a ’’matter” distribution 
outside the galaxy. The existence of the dark radiation term generates an equivalent mass term Mu, which is linearly 
increasing with the distance, and proportional to the baryonic mass of the galaxy, ~ Afs(r/ro). The particles 

moving in this geometry feel the gravitational effects of U , which can be expressed in terms of an equivalent mass. 
Moreover, the limiting value Vt goc of the angular velocity can be obtained as function of the total baryonic mass Mb 
and radius rg of the galaxy as Vt g0 o ~ (2/ \/3)\J GMb/ rg + (l/12\/3 ) (GMs/rg) 3 ^ 2 0. For a galaxy with baryonic 
mass of the order 10 9 Mq and radius of the order of rg ss 70 kpc, we have Vt go a ~ 287 km/s, which is of the same 
order of magnitude as the observed value of the angular velocity of the galactic rotation curves. In the framework of 
this model all the relevant physical parameters (metric tensor components, dark radiation and dark pressure terms) 
can be obtained as functions of the tangential velocity, and hence they can be determined observationally. 

It is the purpose of the present paper to consider vacuum spacetimes on the brane admitting a one parameter 
group of conformal motions, with the vector field generating the motion being non-static and having a more general, 
non-central symmetry. More exactly, we assume that the vector field £ generating the conformal symmetry has the 
general non-static and non-spherically symmetric form £ = £° ( t , r) d/dt+t/ 1 ( t , r) d/dr+£ 2 ( 9 , </>) d/d0+£ 3 {6, 0) d/dcj), 
that is, we also introduce an explicit dependence of £ on the angular coordinates 9 and </>. With this assumption 
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the metric tensor components can be obtained as functions of the conformal factor ip and of the radial coordinate 
r. The gravitational field equations describing the conformally symmetric vacuum brane with general symmetry can 
be reduced to a single differential-integral equation, whose solution can be obtained in terms of the modified Bessel 
functions. The expressions for the metric coefficients, the dark energy and the dark pressure can be expressed in an 
exact parametric form. 

As a physical application of the derived solutions we consider the behavior of a test particle in a stable circular 
orbit on the conformally symmetric brane. Similar to the spherically symmetric non-static case, the angular velocity 
is always a monotonically increasing function of the radial distance and tends, in the limit of large distances, to a 
constant value. This behavior strongly suggest the possibility that a self similar geometry on the brane, generated by 
a general, non-static and non-centrally symmetric group of conformal motions, may provide an explanation for the 
dynamics of the neutral hydrogen clouds at large distances from the galactic center, which is usually explained by 
postulating the existence of the dark matter. 

The present paper is organized as follows. The basic equations describing the spherically symmetric gravitational 
field equations in the vacuum on the static brane are derived in Section II. The general form of the metric tensor 
and the field equations for vacuum brane space-times admitting a one parameter group of conformal motions, with 
non-static and non-central conformal symmetry, are derived in Section III. In Section IV the general solution of the 
gravitational field equations is obtained. The behavior of the angular velocity of a test particle in stable circular 
motion is considered in Section V. We conclude and discuss our results in Section VI. 


II. THE FIELD EQUATIONS FOR A STATIC, SPHERICALLY SYMMETRIC VACUUM BRANE 


We start by considering a 5-dimensional space-time (the bulk), with a single 4-dimensional brane, on which gravity 
is confined. The 4—D brane world g^A is located at a hypersurface (B (J A ) = 0) in the 5—D bulk space-time 

(( 5 ) M, qab), of which coordinates are described by X A , A = 0,1 ,..., 4. The action of the system is given by Jj 


S = 


Sbulk V Sbrane i 


(1) 


where 


and 


Sbulk — 


[ 

Jwm 


! (5) 


2k! 


R -|A 5 ) L m + A 5 


d b X, 


Sbrane 



Sbrane (.9 ot(3 1 Ip') 4“ 'V 


d A x, 


( 2 ) 

(3) 


where = 87 rGs is the 5-D gravitational constant, ^R and ^L m are the 5—D scalar curvature and matter La- 
grangian in the bulk, respectively. x^,g = 0,1,2, 3 are the induced 4—D coordinates of the brane, is the trace 
of the extrinsic curvature on either side of the brane and Lbrane (Sa/Lt^) is the 4—D Lagrangian, which is given by a 
generic functional of the brane metric g a p and of the matter fields ip. In the following capital Latin indices run in the 
range 0,..., 4, while Greek indices take the values 0, ...,3. A 5 and A are the negative vacuum energy densities in the 
bulk and in the brane, respectively. 

The Einstein field equations in the bulk are given by 


( 5 ) G/j = kl^Tu, ( 5 ) T u = -A 5 (5) gu + 6(B) \-X b i5) gu + T u 


(4) 


where 


(5) Tu = -2 


8&L m 

(5( 5 ) g IJ 


(5) 


gu 


(5 )L 


(5) 


is the energy-momentum tensor of bulk matter fields, while T is the energy-momentum tensor localized on the brane 
and which is defined by 


T = —2 

flV — ^ 


8Lbr 


Sg a*" 


g^Lb r 


( 6 ) 


The 8 ( B ) denotes the localization of brane contribution. In the 5-D space-time a brane is a fixed point of the Z 2 
symmetry. The basic equations on the brane are obtained by projection of the variables onto the brane world, because 
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we assume that the gravity on the brane is confined. The induced 4—D metric is gu gu — njnj , where n/ is 
the space-like unit vector field normal to the brane hypersurface G) M. In the following we assume that all the matter 
fields, except gravitation, are confined to the brane. This implies that ^L m = 0. 

Assuming a metric of the form ds 2 = (njnj + gij)dx ! dx J , with njdx 1 = d\ the unit normal to the \ =constant 
hypersurfaces and gu the induced metric on % ^constant hvpersurfaces, the effective four-dinrensional gravitational 
equations on the brane (the Gauss equation), take the form [£j: 


Gfj,v — —hg^ v + k 2 T^ v + k^S^v — E t 


(7) 


where S^ is the local quadratic energy-momentum correction 

- -T«T VOL + ± 9fU , (3 T^T a g - T 2 ) , (8) 

and is the non-local effect from the free bulk gravitational field, the transmitted projection of the bulk Weyl 
tensor Ciajb , Ejj = CiAjBn A n B , with the property Ejj —> E^SjSj as x 0- We have also denoted k 2 = 8ttG, 
with G the usual four-dimensional gravitational constant. 

The four-dimensional cosmological constant, A, and the four-dimensional coupling constant, are given by A = 
k 2 (A ,5 + k 2 A 2 /6) /2 and k 2 = fc|A&/6, respectively. In the limit A,)" 1 —> 0 we recover standard general relativity. 

The Einstein equation in the bulk and the Codazzi equation also imply the conservation of the energy-momentum 
tensor of the matter on the brane, D^T^ 11 = 0, where D v denotes the brane covariant derivative. Moreover, from 
the contracted Bianchi identities on the brane it follows that the projected Weyl tensor should obey the constraint 
D V EJ = ktD v S /. 

The symmetry properties of E^ imply that in general we can decompose it irreducibly with respect to a chosen 
4-velocity field as Q 


E„ v = -fc 4 


U ( 'U'n'U'i/ + ^ ) + PU V + 2 Q^u v) 


(9) 


where k = k$/k 4 , h^ = g^ + u^u v projects orthogonal to the ’’dark radiation” term U = —k i E lxl/ u^u 1 ' is a 
scalar, Q M = k A h^E a p a spatial vector and = —fc 4 E a p a spatial, symmetric and trace-free 

tensor. 

In the case of the vacuum state we have p = p = 0, = 0 and consequently = 0. Therefore the field equations 

describing a static brane take the form 


Riiu — —E^v + 4A, (10) 

with the trace R of the Ricci tensor R satisfying the condition R = Rfj, = E£ = 0. 

In the vacuum case E^ satisfies the constraint D^E^ = 0. In an inertial frame at any point on the brane we have 
u 11 = Sq and h^ =diag(0,1,1,1). In a static vacuum Q^ = 0 and the constraint for E^ takes the form EH 

^D li U+^UA li + D v P ltv + A v P ltv = 0, (11) 

where D M is the projection (orthogonal to ip) of the covariant derivative and A M is the 4-acceleration. 

In the static spherically symmetric case we may chose A M = A(r)r^ and P M „ = P(r) (r^ry, — |/i M1 ,), where A(r ) and 
P(r) (the ’’dark pressure”) are some scalar functions of the radial distance r, and r^ is a unit radial vector jlCj . 

We chose the static spherically symmetric metric on the brane in the standard form 

ds 2 = -e^ r) dt 2 + e A(r W 2 + r 2 (. id 2 + sin 2 Odcf 2 ) . ( 12 ) 

Then the gravitational field equations and the effective energy-momentum tensor conservation equation in the 
vacuum take the form 0,0 


—e 


-A 




48t tG 
/c 4 Xb 


U + A, 


(13) 


e 


-A 



1 

^2 


167rG 

k 4 X b 


(U + 2 P) - A, 


r 


(14) 
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e-M*/' + T + 


j/2 v '_\' j/\'\ 327 tG 




(15) 


U' + 2P' 


6 P 


2U + P r (2U + P) ’ 

where we denoted ' = d/dr. In the following we denote (3 = 167rG/fc 4 Af,. 


(16) 


III. FIELD EQUATIONS ON A VACUUM BRANE WITH NON-STATIC AND NON-CENTRAL 

CONFORMAL SYMMETRY 


The system of the field equations for the vacuum on the brane is under-determined. A functional relation between 
the dark energy U and the dark pressure P must be specified in order to solve the equations. An alternative method, 
which avoids ad hoc specifications, is to assume that the brane is mapped conformally onto itself along the direction 
£, so that 


— 9^v, a£ + + 9n — VdUj' i 


(17) 


where 0 is the conformal factor. In the following we shall assume that 0 is a function of the radial coordinate only, 
0 = 0(r). As for the choice of £, in the framework of the standard general relativity, Herrera and Ponce de Leon 
assumed that the vector field generating the conformal symmetry is static and spherically symmetric, 


Using this form of the conformal vector in Eqs. m one obtains 

e = A ,e = ^re X p(~ 


and 


0 (r) = B exp ( — ^ , exp (i/) = C 2 r 2 , 


(18) 


(19) 


( 20 ) 


respectively, where A , P, C are constants. A may be set to zero since Ad/dt is a Killing vector and B may be set to 
1 by a rescaling £ —> P -1 £, 0 —> P _1 0, which leaves Eqs. G3 invariant. 

This form of £ gives the most general £ invariant under the Killing symmetries, that is [<9/3i,£] = 0 = [X a ,£], 
where X Q generates SO (3). The corresponding form of the metric, obtained by imposing static conformal symmetry, 
has been used to study the properties of compact anisotropic general relativistic objects Q and to investigate the 
properties of charged strange stars ,15]. The general solution of the vacuum brane gravitational field equations for 
this choice of £ has been obtained in 1]. 

A more general conformal symmetry has been proposed by Maartens and Maharajah 0, which generalizes the 
isotropic conformal vector td/dt + rd/dr of the Minkowski space-time, but weakens the static symmetry of £ in Eq. 

m- 




( 21 ) 


With this choice the conformal vector is non-static and spherically symmetric. With the assumption that the 
conformal factor 0 is static, 0 = 0 (r) and with this form of £ Eqs. ® give 0: 


£° = A + 


1 k 
2B* 


( 22 ) 


= C 2 r 2 exp 




(23) 
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ip = Be~ x/2 , 


(24) 


where k is a separation constant and A, B and C are integration constants. Without any loss of generality we can 
chose A = 0. Thus in this model we obtain for the vector field £ the expression 


1 k ^ d rip (r) d 

2 B l di + 2 dr' 


(25) 


In the framework of the standard general relativity the static solutions of the gravitational field equations with 
metric tensor components given by Eq s. (1231) and <EH) have been obtained, for anisotropic and charged compact 
objects, by Maartens and Maharaj [Hj. All the solutions of the vacuum gravitational field equations on the brane 
having the conformal symmetry generated by the vector £ given by Eq. ( 1231 ) have been obtained in @|. 

In the following we consider a generalization of the conformal vector £ by relaxing the condition of the spherical 
symmetry. Hence we assume a more general form of the vector field £, generating the conformal symmetry, namely 

{=<<, e |+c (*, e |:+1 2 («■« ^+ e («,« (20 


We still assume that the conformal factor is spherically symmetric, that is, we take ip 
form of £, Eqs. (THT) give 


ip (r). With the use of this 


^+ 2 ^ = (27) 

X'e+^=iP( r ), (28) 


g 1 . d£ 2 _ ^(r) 

r dO 2 


(29) 


From Eq. i!231i it follows that 


? , COS0 2 ^ i>(r) 

r sin 0 dtp 2 


g 1 (r) 


r (ip — ex) 
2 


(30) 


(31) 


g 2 (M) = |* + 


d<p 


(32) 


where a is an arbitrary separation constant and F an arbirary function of the azimuthal angle <p. Then Eq. pl gives 


£ 3 (0, (p) = | (1-0 cot 0)<p -cot OF (</>) + G(0), (33) 

where G ( 0 ) is an arbitrary integration function. From Eq. 1271) we obtain 

£°(t) = |t + A, (34) 

where A: is a separation constant and A is an integration constant, which can be taken to be zero without any loss of 
generality, A = 0. Hence the metric tensor component exp ( v ) is given by 


exp (y) = C 2 r 2 exp 




dr 

r (ip — a) \ 


(35) 
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where C is an arbitrary integration constant, 
form 

exp (A) = 


Finally, the metric tensor component exp (A) can be obtained in the 


B 2 


O - «)' 


■ exp 


2 a 


dr 


r {ip — a) 


(36) 


where B is an arbitrary constant of integration. 

Therefore the requirement of the conformal symmetry on the brane completely fixes the form of the metric tensor, 
as a function of the conformal factor ip. As for the vector field £, it is given by 


„ k d 

t = ?m 


r {ip — a) d 
2 dr 


a . 


dF{<P) 


dp 


d_ 

dd 


a 

12 


(1 — 0 cot 6)<p — cot 9F {(p) + G (0) 


d_ 

dp 


(37) 


With the use of the representations given by Eqs. (l35l) and for the metric tensor components, the gravitational 
field equations CS-C0 on the vacuum brane with general, non-central conformal symmetry, take the form 


{4> ~ a) 2 

B 2 


{ip — a) 2af ^^ 
B 2 


(iP - a ) 2 2a f 
B 2 


2ip' 


= 2q / TlP-a) 



-3 

s 

r 2 {ip 

dr 

r{ip — Q 

T [ 2 {a - k) 

3 

-1- 


1 

to 

A 

' r 2 

2 ip' 

2k 



2a 1 

7 d-o 


-|-y — 3 j3U + A, 


- o = P{U + 2P) — A, 


k{a — k ) 1 

r(ip-a) r 2 {ip - a) r 2 _ a f r 2 


= (3 (U — P) — A. 


(38) 


(39) 


(40) 


IV. GENERAL SOLUTION OF THE GRAVITATIONAL FIELD EQUATIONS FOR THE VACUUM 
BRANE WITH NON-CENTRAL CONFORMAL SYMMETRY 

By multiplying Eq. m with 2, adding the obtained equation to Eq. cm and equating the resulting equation with 
Eq. cm gives the following differential-integral equation, satisfied by the conformal factor ip (r): 

3 r {ip — a) ip' + 3 ip 2 — 3 {k + 2a) ip + 3a 2 + 2 ak + k 2 — B 2 e 2a ^ _ 4A B 2 e 2a ^ r w- Q > r 2 = 0. (41) 

In order to solve Eq. m we introduce a new variable u defined as 

dr 


lnw = a 


r {ip — a) 


(42) 


In terms of u Eq. <ED can be expressed in the form 

—3— u 2 ^~o +3 a {a — k)^-— + k 2 — ka — B 2 u 2 — 4A B 2 u 2 r 2 = 0. (43) 

r u 6 r u 

With the use of the identities 1/u' = dr/du and u"/u' 3 = —d 2 r/du 2 , Eq. 14311 can be transformed into the following 
differential equation: 

u2 ~T~5 ~ mt *T— U 2 u 2 — n 2 ) r — 12lAu 2 r 3 = 0, (44) 

du 2 du 

where we denoted 


to = 



l 2 



k {k — a) 
3a 2 


(45) 


In the following we assume that we can neglect the nonlinear term in Eq. m . This approximation is valid in the 
case of very small numerical values of the parameters l and A and for small values of r. Therefore the basic equation 



































describing the geometry of the vacuum brane with geometric self-similarity with the effect of the cosmological constant 
ignored is given by the following linear, Bessel type differential equation: 


u2 ~^~2 ~ mu ~j — ( l 2 u 2 — n 2 ) r = 0 , 


dr 

du 


The general solution of Eq. cm is given by 

r (u) = u s [Clip ( lu ) + C 2 K p (Zu)] , 


(46) 


(47) 


where C\ and C 2 are arbitrary constants of integration, s = (l + m)/2 = k/2a , p = yj (1 + m) 2 — 4n 2 /2 = 


\Js (2 — s)/3, and I p ( x) and K p (x ) are the modified Bessel functions of the first and second kind, respectively, 
satisfying the differential equation x 2 y" + xy' — (x 2 + p 2 ) y = 0 E2- For all values of p , I p [x) and K p (x) are linearly 
independent. I p (x ) and K p (x) are real and positive when p > —1 and x > 0. Since by definition l > 0, it follows that 
the parameter u must be non-negative, u > 0. On the other hand the condition that p must be a real number, p £ R, 
imposes the constraint 0 < s < 2 on the numerical parameter s. 

To obtain a correct physical description, the function r(u) must be an increasing function of the variable u and 
obey the conditions lim u _>o r(u) = 0 and hin^^oo r(u) = oo, respectively. For small values of the arguments and for 
fixed p the modified Bessel functions behave like I p (x) ~ ( x/2) p /T(p + 1) and K p {x) ~ T (p) ( x/2)~ p /2, respectively, 

G3- 


’ t z x e z dt 171. Hence in the limit of small u we obtain 


where T is the Euler gamma function, defined as T (z) = f 0 
lim„_,o r(u) ~ u s C\ (l/2) p u p /T(jp + 1) + C 2 V (p) (1/2)~ p u~ p /2 . In order to obtain lim„^o r(u) = 0, the condition 

s > p must necessarily hold. This condition is satisfied for values of s so that s > 1/2. 

In the opposite limit of large values of the argument, the modified Bessel functions I P (x) and K p (x) have the 
asymptotic behaviors I p (x) ~ exp (x) /\/2itx and K p {x) ~ tJtt/2x exp (— x), respectively ^3- Hence lim^^oo r(u) = 
lim. 


(Ci/ 


2ttIju s 1 / 2 exp(lu) 


oo. To obtain the correct form of the limit it is also necessary that G\ > 0. 

The numerical values of the parameter s are restricted to the range s £ (1/2, 2]. 

Therefore the general solution of the gravitational field equations on the vacuum brane with non-static and non¬ 
central conformal symmetry can be expressed in an exact parametric form, with u taken as parameter, and given by 


e" = CV (1 ~ s) [Cil p (lu) + C 2 K p ( lu)f 


(48) 


A 3l 2 u 2s [C 1 I p (lu)+C 2 K p (lu)Y 

6 ~ 7dir 

V du ) 


(49) 


^ u =M l+ h 


d , V 2(2a-l) ( d , 

— In r H-— In r 

du J u \du 


2 n 2 


(50) 


3 /3P = 


3l 2 r 2 


4 ( — In r | — 
du 


2 ( 2 s- 1 ) ( d 


du 


lnr ) + l" — 


r-2 ‘ 


(51) 


The variation of the metric coefficient exp(i/) jC 2 C\ is represented, as a function of r/C 2 , in Fig. 1 

For the chosen values of the parameters exp(v) is a decreasing function of the radial distance r. The variation of 
the metric function exp(A)/31 2 is represented in Fig. 2. 


The metric function exp(A) is a monotonically increasing function of the coordinate r. In the limit of large values 
of r, exp(A) tends to a constant value. The variation of the dark radiation term 3 PC 2 U is represented, as a function 
of r/C 2 , in Fig. 3. 
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FIG. 1: Variation, as a function of the radial distance r/C 2 ^ of the metric coefficient exp(^) jC 2 C\ (in a logarithmic scale) for 
a static, conformally symmetric vacuum space-time on the brane, for l = 10~ 9 , Ro = C 1 /C 2 = 10~ 5 and different values of s: 
s = 1.1 (solid curve), s = 1.2 (dotted curve), s = 1.3 (dashed curve) and s = 1.4 (long dashed curve). 



r/C 2 

FIG. 2: Variation, as a function of the radial distance r/C?, of the metric coefficient expA/31 2 (in a logarithmic scale) for a 
static, conformally symmetric vacuum space-time on the brane, for l = 1CF 9 , Ro = C 1 /C 2 = 1CP 5 and different values of s: 
s = 1.1 (solid curve), s = 1.2 (dotted curve), s = 1.3 (dashed curve) and s = 1.4 (long dashed curve). 


The dark radiation term is positive for all values of the radial coordinate r, U(r) > 0, Vr £ (0, 00 ). In the limit of 
large r, U tends to zero, lim^oo U(r) = 0. The variation of the absolute value of the dark pressure 3/3C| |P| as a 
function of r/C 2 is represented in Fig. 4. 


In the present model the dark pressure is negative, satisfying the condition P(r) < 0, Vr £ (0,oo). In the large 
time limit, similar to the dark radiation term, the dark pressure also tends to zero, linv^oo P(r) = 0. 

In the limit of small r, corresponding to small values of the parameter u, we have r « (C 2 / 2 )T(p) (l/2)~ p u s ~ p , 
giving 


u 


2 

C^p) 



(52) 


Therefore in the limit of small radial distances the metric coefficients e" and e x behave as 

2(l-2s) 

s ~ p 2(1 —s—p) 

r s -p , r —> 0 , 


C 2 


C2T(p) V 2 


(53) 


and 


3 1 2 


c 2 r(p) 


1 + 


r( P +i) 

TTpT 


(l \ 2 r(p-i) 


2 

c 2 r( P ) 


(k ) 1 


(54) 
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FIG. 3: Variation, as a function of the radial distance r/C‘ 2 , of the dark radiation term 3/3(711/ (in a logarithmic scale) for a 
static, conformally symmetric vacuum space-time on the brane, for l = 1CV 9 , Ro = C\/C 2 = 10~ 3 and different values of s: 
s = 1.1 (solid curve), s = 1.2 (dotted curve), s = 1.3 (dashed curve) and s = 1.4 (long dashed curve). 



FIG. 4: Variation, as a function of the radial distance r/C 2 , of the dark pressure term 3/3(71 13 d ! for a static, conformally 
symmetric vacuum space-time on the brane, for l = 10 9 , Ro = C 1 /C 2 = 1CV 3 and different values of s: s — 1.1 (solid curve), 
s = 1.2 (dotted curve), s = 1.3 (dashed curve) and s = 1.4 (long dashed curve). 


respectively. 

In the limit of large r the dependence of the metric functions on the radial distance can be given only in a parametric 
form. Since for u —> 00 we have r = (^Ci/y/2nCj u s-1 / 2 exp(lu), the asymptotic behavior of exp (v) and exp (A) is 
given, for large r, in the parametric form, with u taken as parameter 


C 


y/2nl 


1 u s ~ 1 ^ 2 e lu , e" 


^ '-"'I 1— 2<t 0.111 

r u 2 Se 2 i v->°°’ 


and 


C 1 


■.u s ~ 1 / 2 e lu .e x 


2„.2 


3 ru 


00 . 


\p2jd (s + luf 

For r —> 0 the radial distance dependence of the dark radiation U is given by 


(55) 


(56) 


where 


3/317 (r) 


1 

3?’ 2 



Ku 


2 n 2 + [s + r(p + l)/r(p )] 2 + 2(2s - 1) [s + F(p + 1)/T(p)} 


r- 


2 (i\P 
c 2 r( P ) ( 2 ) 


2 

s-p 


(57) 


(58) 
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In the limit of large r, r —> oo, the dark radiation behaves as 


3 PU{r) « ^,r -> oo. 


In the limit of small r we find for the dark pressure 


where 


K P = 


3/3P(r) « —f 1 + —) , r —> 0, 


n 2 - 4 [s + r(p + i)/r ( P )] 2 + 2(2 s -1) [s + r(p + i)/r(p)] 


6/ 2 


2 ( l\P 

c 2 r( P ) 




For r —> oo the dark pressure has the asymptotic limit 


1 


3 /3P(r) « , r oo. 

Therefore at infinity the dark pressure and the dark radiation obey the equation of state 

U + 2P = 0. 


(59) 

(60) 

(61) 

(62) 

(63) 


V. STABLE CIRCULAR ORBITS IN VACUUM BRANE SPACE-TIMES WITH GENERAL 

CONFORMAL SYMMETRY 


As a physical application of the conformally symmetric brane metric m generated by a vector field £ with arbitrary 
symmetry we consider now the problem of constructing stable circular timelike geodesic orbits. The motion of a test 
particle in the gravitational field can be described by the Lagrangian [3 




2 idfiy 
r '.*) • 


(64) 


where we denoted by r the affine parameter along the geodesics. In the timelike case r corresponds to the proper 
time. The equation giving the tangential velocity of the body has been derived in 0 and j]j|, and is given by 


r v 


(65) 


Thus, the rotational velocity of the test body is determined by the metric coefficient exp (v) only. 

In the case of the motion of a test particle in a conformally symmetric, static spherically symmetric space-time, 
with a general non-static vector field generating the symmetry, the metric coefficient exp ( v ) is given by Eq. lt23l) . 
Therefore as a function of the conformal factor (/> the angular velocity Vt g is given by the simple expression 


v t 9 = 1 ^ 


*(2s- 1) 
— a 


( 66 ) 


Eq. Ki6l> gives a simple physical interpretation of the conformal factor ift in terms of the tangential velocity, 
if) = a 1 + (2s — 1) (l — Vf g ) 1 . Hence the conformal factor is proportional to the tangential velocity of a test 

particle in the gravitational field of a brane with general conformal symmetry. As a function of the parameter u, vt g 
is given by 


V tg = 1 ^ 


2 s — 1 

u £ h ^ r 


(67) 


Together with Eq. (IT71) . Eq. m gives the parametric representation of the tangential velocity of a test particle 
in a stable circular orbit on the radial distance r. By using the parametric dependence of r we obtain the following 
explicit representation for vt g as a function of u: 


Vt g (u) = l- 


2 s- 1 


„ I lu Ci[Ip-i(lu)-\-Ip+i(lu)]-\-C 2 [Kp — i(lu)-\-K p -\-i(lu)\ 
CiI p (lu)+C 2 K p (lu ) 


( 68 ) 


2 
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FIG. 5: Variation, as a function of the parameter r/C 2 , of the tangential velocity vt g of a test particle in a stable circular orbit 
in a conformally symmetric vacuum space-time with a general conformal vector on the brane, for l = 10~ 3 ' 2 , Ro = C 1 /C 2 = 10 s 
and different values of s: s = 1.1 (solid curve), s = 1.101 (dotted curve), s = 1.102 (dashed curve) and s = 1.103 (long dashed 
curve). 


The variation of Vt g as a function of the radial distance r is represented, for different values of the constant s, in 
Fig. 5. 


The tangential velocity of a particle in a stable circular orbit in a conformally symmetric geometry on the brane, 
with arbitrary conformal factor, is a monotonically increasing function of the distance. 

In the limit of small radial distances, by using the approximate expressions of the Bessel functions for small values 
of the argument we obtain the following explicit dependence of the tangential velocity on the radial distance r: 

<i, if) ~ 1----- ~ —- ------, r -► 0. (69) 

s + T? + (&) * * (s) S " P T i p “ (P) 


The angular velocity is a monotonically increasing function of r. For small values of r, one can neglect the term con¬ 
taining the radial coordinate with respect to the other terms, giving v\ g sa 1 — (2s — 1)/ [s + T(p + 1)/T (p)) =constant. 
Therefore at small distances from the origin the angular velocity of a test particle moving on the conformally symmetric 
brane has a constant value. 

In the limit of large distances we obtain the following parametric representation for vt g . 


r 




1 - 


2 s — 1 
s + lu 


(70) 


In the limit of large r the tangential velocity of a test particle in a stable circular orbit in a vacuum brane with 
general conformal symmetry tends to the speed of light c. This type of behavior, corresponding to an increase in the 
velocity of test particles in a stable circular orbit, is very similar to the behavior of hydrogen clouds outside spiral 
galaxies, and which is usually attributed to the presence of the dark matter. Therefore, based on this analogy, we 
shall assume that the conformally symmetric solution of the gravitational held equations on the brane describes the 
gravitational held outside a galaxy. 

The metric coefficient exp(A) can be expressed as a function of the tangential velocity and the parameter u as 


e 


A 




(71) 


The held equation m can be immediately integrated to give the following representation for the metric tensor 
component exp (—A), 


where 


c -x _ 1 IGMuir) . 


Mij ( r ) = / u i r ) r2dr - 


(72) 


0 


(73) 
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is the mass corresponding to the dark radiation term U(r) (the ’’dark mass”). By using Eq. 1491) we obtain the 
following parametric representation for Mjj(r): 


Mu(r) = 


2 G 


1 - 


1 ( d 


3Z 2 \ du 


■In; 


(74) 


The dark mass can be also expressed in terms of the tangential velocity of a particle in a stable circular orbit as 


M„(r) = ^ 


t (2» ~ if _1 


u 2 (l — V? 


tg) 


(75) 


In the limit of small r the dark mass behaves like 


Mu ( r ) ~ 2G i 1 3 1 2 


s + 


r (p + 1 ) 
r (p) 


-i 2 r 


C 2 T(p) \2 


pi 


-p i 


, r —0. 


(76) 


In the limit of large r (u —> oo) we obtain 


Mtr(r) = 3 a 


(77) 


Therefore for large r the dark mass is linearly increasing with the distance to the galactic center. 


VI. DISCUSSIONS AND FINAL REMARKS 

In the present paper we have obtained a class of conformally symmetric solutions of the vacuum field equations 
in the brane world model, under the assumption of a non-static and non-central conformal symmetry, and we have 
discussed some of their physical properties. To obtain the solution we have made the crucial assumption of ignoring the 
possible effect of the cosmological constant on the geometrical structure of the vacuum on the brane. Mathematically, 
this means that the model can correctly describe the dynamics of particles at distances r from the galactic center 
satisfying the condition 


r « 


l 2 u 2 - n 2 
\2lku 2 


(78) 


Taking into account the smallness of the value of the cosmological constant, the present model can be safely applied 
to describe the motion of particles outside galaxies. On the other hand, in the limit of large r, both the dark radiation 
U and the dark pressure P are decreasing functions of r, and at infinity the geometry and the dynamical behavior of 
the particles is determined by the cosmological constant only. 

The behavior of the metric coefficients and of the angular velocity in the solution we have obtained depends on four 
arbitrary constants of integration Z, s, C\ and C 2 . Their numerical values can be obtained by assuming the continuity 
of the metric coefficients exp (A) and exp [y) across the vacuum boundary of the galaxy. For simplicity we assume 
that inside the ’’normal” (baryonic) luminous matter, with density ps, which form a galaxy, the non-local effects of 
the Weyl tensor can be neglected. We define the vacuum boundary tq of the galaxy (which, for simplicity, is assumed 
to have spherical symmetry) by the condition pb (ro) ~ 0. 

Therefore at the vacuum boundary the metric coefficients are given by exp (y) = exp (—A) = 1 — 2GMb/tq, where 
Mb =47r f’° ps (r) r 2 dr is the total baryonic mass inside the radius ro- For simplicity we shall also assume that the 
tangential velocity at the galactic boundary can be approximated by its Newtonian expression, v 2 g (ro) ~ GMb/tq. 

The continuity of the radius r through the surface r = tq gives the first condition which must be satisfied by the 
integration constants: 


r 0 = u s 0 [Clip ( lu 0 ) + C 2 K p ( lu 0 )] , (79) 

where uq is the value of the parameter u corresponding to the vacuum boundary of the galaxy r = tq. The continuity 
of exp (y) at r = ro gives 


C 2 r 2 0 u 0 {1 - 2s) = l-2GM B /r 0 . 


(80) 
















14 


The continuity of the tangential velocity of a test particle across the galactic boundary fixes the value of the 
derivative dr/du for u = Uq\ 


dr\ ro 2s — 1 

du )n=uo 


Then the continuity of exp (A) gives 

4-^) (1 - 2GM B /r 0 ) = i 


2 s- 1 

l 


(81) 


(82) 


Thus the continuity of the metric potentials gives the following compatibility condition, which must be satisfied by 
the constants C, l and s: 


C 2 


1 2s — 1 


2(1 —2s) 


3 l 


1_ f _ GM b \ 2{1 ~ 2s) f _ 2 GM B \ l ~ 2s 

r ol r 0 ) V r 0 ) 


(83) 


As a possible physical application of the obtained solutions we have considered the behavior of the angular velocity 
of a test particle in a stable circular orbit on the brane with non-static and non-central symmetry. The conformal 
factor ip, together with two constants of integration, uniquely determines the rotational velocity of the particle. The 
angular velocity is always an increasing function of the radial distance r, which from a physical point of view can 
be considered as the distance from the galactic center. This behavior is independent on the numerical values of the 
parameters (separation and integration constants) of the model. In the limit of large radial distances the angular 
velocity tends to a constant value, which in the present case is the speed of light c. This general behavior is typical for 
massive particles (hydrogen clouds) outside galaxies. Thus the rotational galactic curves can be naturally explained 
in brane world models, without using the concept of dark matter. 

It has long been known that Newtonian or general relativistic mechanics applied to the visible matter in galaxies 
and clusters does not correctly describe the dynamics of those systems. The rotation curves of spiral galaxies |13j are 
one of the best evidences showing the problems Newtonian mechanics and/or standard general relativity has to face 
on the galactic/intergalactic scale. In these galaxies, neutral hydrogen clouds are observed at large distances from 
the center, much beyond the extent of the luminous matter. Assuming a non-relativistic Doppler effect and emission 
from stable circular orbits in a Newtonian gravitational field, the frequency shifts in the 21 cm line hydrogen emission 
lines allows the measurement of the velocity of the clouds. 

Observations show that the rotational velocities increase near the center of the galaxy and then remain nearly 
constant at a value of v tg0 o ~ 200 km/s [Til . This leads to a mass profile M(r ) = rvf goo /G. Consequently, the 
mass within a distance r from the center of the galaxy increases linearly with r, even at large distances where very 
little luminous matter can be detected. This behavior of the galactic rotation curves is explained by postulating the 
existence of some dark (invisible) matter, distributed in a spherical halo around the galaxies. The dark matter is 
assumed to be a cold, pressureless medium. There are many possible candidates for dark matter, the most popular 
ones being the weekly interacting massive particles (WIMP). 

However, despite more than 20 years of intense experimental and observational effort, up to now no non-gravitational 
evidence for dark matter has ever been found: no direct evidence of it and no annihilation radiation from it. Moreover, 
accelerator and reactor experiments do not support the physics (beyond the standard model) on which the dark matter 
hypothesis is based. 

Therefore, it seems that the possibility that Einstein’s (and the Newtonian) gravity breaks down at the scale of 
galaxies cannot be excluded a priori. Several theoretical models, based on a modification of Newton’s law or of general 
relativity, have been proposed to explain the behavior of the galactic rotation curves (2jj. In the framework of the 
brane models, the role of the dark matter is played by the dark mass Mjj{r ), which is the mass associated to the dark 
radiation term, having its physical origin in the five-dimensional bulk. As can be seen from Eq. GZD, the dark mass 
is indeed linearly increasing with the distance to the galactic center. Thus the dark mass has a similar behavior to 
that observed in the case of the ’’dark matter” around the galaxies. 

The main advantage of the brane world models as compared to the alternative explanations of the galactic rotation 
curves and of the dark matter is that it can give a systematic and coherent description of the Universe from galactic 
to cosmological scales. On the other hand, in the present model all the relevant physical quantities, including the 
dark energy and the dark pressure, which describe the non-local effects due to the gravitational field of the bulk, are 
expressed in terms of observable parameters (the baryonic mass and the radius of the galaxy). Therefore this opens 
the possibility of the testing of the brane world models by using astronomical and astrophysical observations at the 
galactic scale. A systematic comparison between the predictions of the model discussed in the present paper and the 
observational results will be considered in detail in a future publication. 
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